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We discuss the production of W bosons in association with three jets at the LHC. We investigate 
how next-to-leading order QCD corrections modify basic kinematic distributions of jets and leptons. 
We also address the magnitude of NLO QCD effects in -I- 3 jet observables, relevant for SUSY 
searches at the LHC. 



I. INTRODUCTION 

A good understanding of complicated multi-particle 
processes is important for the LHC physics. To achieve 
this goal it is useful to have next-to-leading order (NLO) 
QCD predictions for such processes (see e.g. Ref. In 
the past, three final state particles was the highest multi- 
plicity for which NLO QCD computations were feasible, 
but this changed this year when four groups [1, H, 0, H, 
SSIl] reported first results on NLO QCD corrections to 
processes with Jour particles in the final state. To arrive 
at these results a variety of methods including highly- 
refined Passarino- Veltman reduction algorithm , 
Ossola-Pittau-Papadopqulos (OPP) method (Til, [l^ and 
unitarity techniques [H, [13, IH, [3, [Hi were used. 
Successful completion of these computations and a large 
number of one-loop amplitudes with six and more exter- 
nal particles computed recently ESUHil,!!!,!!!!!!], 
provides a proof of principle that reliable description of 
many 2 — > 4 processes is now within reach. 

A major reason for extending leading order results to 
next-to-leading order is a significant reduction of the un- 
physical dependence on factorization and renormaliza- 
tion scales in NLO QCD cross sections and distributions. 
Such reduction is very important, especially for high mul- 
tiplicity processes, where the unphysical dependence on 
scales can be significant. Indeed, for such processes, the 
scale dependence is amplified by the high power of the 
strong coupling constant. For a cross-section involving n 
jets iT„ ~ a"(/^), so that small changes in the renormal- 
ization scale /i lead to large changes in the corresponding 
cross sections 

-2n/3oa,a„, /3o = (llA^c - 2n/)/(127r). (1) 
There are many cases where NLO computations reduce 



the dependence of cross sections on unphysical scales to 
10 — 20% which, given typical cross section uncertainties 
at leading order (LO) of about 50%, is a great success. 
However, even if such scale-independence is observed for 
the total-cross section, it is not always possible to claim 
that a particular process is described reliably, for the 
purpose of LHC phenomenology. The reason is the dual 
role that many complicated processes at the LHC will 
play. Indeed, depending on the cuts on the final states, 
processes that involve e.g. top quarks and/or electroweak 
gauge bosons and QCD jets are treated as either primary 
signals or unwanted backgrounds and very different cuts 
are applied to final states in the two cases. We will re- 
fer generically to these cuts as "signal" or "background" 
cuts, respectively ^ These cuts force final state particles 
to live in different regions of phase-space, so that a'priori 
it is not possible to relate QCD corrections to the same 
process subject to either "signal" or "background" cuts. 

A candidate procedure for dealing with LHC processes 
which have this dual role is as follows. One starts the 
study of these processes by applying signal cuts and us- 
ing the resultant data set to refine the theoretical tools, 
e.g. Monte Carlo event generators. Once a good under- 
standing of a given process is achieved with signal cuts, 
one uses the refined tools to extrapolate to a different 
kinematic situation, specified by the background cuts. 
Unfortunately, the reliability of this extrapolation is not 
assured. The purpose of applying the background cuts 



Since we are mainly concerned with the process W + 3 jets, we 
call "signal cuts" the ones where this process is considered a 
signal, and "background cuts" the ones for which this process 
is an unwanted background to some other New Physics signal. 
We caution the reader that these terms are often used in the 
literature in exactly the opposite way. 
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FIG. 1: The dependence of the + 3 jet inclusive produc- 
tion cross section at the LHC on the factorization and renor- 
malization scale /i. All cuts and parameters are described in 
the text. The leading color adjustment procedure is applied. 
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FIG. 2: The dependence of the W~ + 3 jet inclusive produc- 
tion cross section at the LHC on the factorization and renor- 
malization scale fi. All cuts and parameters are described in 
the text. The leading color adjustment procedure is applied. 



is to suppress, as far as possible, the very kinematic con- 
figurations that are allowed by signal cuts. Therefore, 
such an extrapolation can only work if the influence of 
kinematics on QCD radiative effects is correctly captured 
by the available tools. Since only relatively simple theo- 
retical tools, such as leading order parton integrators or 
parton showers, are currently available for complicated 
final states, the modeling of the radiative effects is only 
approximate. On the other hand, if a NLO QCD com- 
putation is available, such an extrapolation can be done 
with a smaller ambiguity since all the relevant scales are 
generated dynamically in NLO computations, largely in- 
dependent of the choices made initially. For cases with 
complicated kinematics, this is clearly indispensable. 

In this paper we discuss and illustrate this issue, taking 
the production of W bosons in association with three jets 
at the LHC as an example. For definiteness, we choose 
to consider proton- proton coUisions at ^/s = 10 TeV [26| . 
However, we do not aim to describe -t- 3 jet produc- 
tion at the LHC in all possible detail, since knowledge 
of the exact experimental setup would be required. In- 
stead, we look for and try to understand differences be- 
tween NLO and LO QCD results for basic observables, 
for the case of signal cuts. We point out that it is not al- 
ways clear which leading order predictions should be used 
in those comparisons since different choices of renormal- 
ization and factorization scales affect the leading order 
predictions strongly. We therefore compare our results 
to a variety of leading order predictions including most 
advanced ones, where matrix element computations are 
matched to parton showers. 

We note that NLO QCD corrections to W+3 jet pro- 
duction at the LHC have been studied in great detail 
recently in Ref. @, mostly for signal cuts. We have 
checked a number of results for production cross- 
sections at the LHC, reported in that reference, and 
found agreement within a few percent in all cases con- 



sidered. These small differences are compatible with the 
fact that our calculation employs the leading color ap- 
proximation with color adjustment procedure, explained 
in detail in Ref. 5], whereas computation in Ref. ac- 
counts for complete color dependence. 

We also discuss QCD corrections to background cuts, 
studied by the ATLAS [13, M, HI] and CMS collabo- 
rations [20, HH for SUSY searches at the LHC. Such 
analyses often assume that Standard Model backgrounds 
can be measured in SUSY-free regions to fix normal- 
izations and then employ LO computations to extrap- 
olate to kinematic regions where supersymmetric signal 
is expected. Hence, an implicit assumption in those 
analyses is that LO distributions have correct shapes 
and that higher-order QCD effects provide a kinematic- 
independent renormalization. We are now in position 
to check these assumptions with the explicit NLO QCD 
computation of W^-|-3 jet process for typical ATLAS and 
CMS cuts. 

The remainder of the paper is organized as follows. In 
Sectionini we discuss M^-l-B jet production for signal cuts 
at the LHC. In Section |TTT] we study 14^-1-3 jet production 
as a background to SUSY searches for two typical sets of 
cuts close to those suggested by the ATLAS and CMS 
collaborations. In Section HVl we present our conclusions. 



II. STUDY OF W + 3 JET PROCESS 

In this Section, we discuss NLO QCD effects in lU -f 3 
jet production for a set of cuts, designed to study the W 
production in association with jets. We follow Ref. Q 
closely and perform calculations in the leading color ap- 
proximation. The calculation relies heavily on the frame- 
work provided by MCFM [33| and uses one-loop ampli- 
tudes computed in l23. We employ the Catani-Seymour 
dipole subtraction [33| to compute real emission correc- 
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FIG. 3: The transverse momentum distribution of the leading 
jet for + 3 jet inclusive production cross section at the 
LHC. All cuts and parameters are described in the text. The 
leading color adjustment procedure is applied. 



FIG. 4: The transverse momentum distribution of the lead- 
ing jet for + 3 jet inclusive production cross section at 
the LHC. All cuts and parameters are described in the text. 
The leading color adjustment procedure is applied. All LO 
distributions are rescaled by constant factor, to ensure that 
the LO and NLO normalizations coincide. 



tion; details of our implementation are given in We 
use the leading color adjustment procedure described in 
that paper to correct for deficiencies of the leading color 
approximation, to the extent possible. We note that pro- 
duction cross-sections for W'^ and W~ at the LHC are 
not the same; we have chosen to discuss the case of 
production almost everywhere in this paper. We do, how- 
ever, show results for the W~ + 3 jet production cross- 
section at the LHC in dependence of factorization and 
renormalization scales. 

We begin by summarizing all the relevant cuts and 
input parameters that are employed in the computation. 
We take the LHC center-of-mass energy to be 10 TeV. We 
require that the transverse momentum and pseudorapid- 
ity of the three jets satisfy ptj > 30 GeV and \T]j\ < 3. 
We consider the leptonic decay of the W to electron (or 
muon) and employ the following restrictions on lepton 
transverse momentum, missing transverse energy, lepton 
rapidity and W^-boson transverse mass, pT,e > 20 GeV, 
> 15 GeV, |77e| < 2.4, > 30 GeV. We 

do not apply an isolation cut on the leptons. To de- 
fine jets, we use the SISCone jet-algorithm [U with 
R = a/Att^ + Ai^ = 0.5 and merging parameter / = 0.5. 

We consider the production of on-shell bosons, 
that decay into a pair of massless leptons. Finite width 
effects are about 1%; they tend to decrease the cross 
section. The CKM matrix is set equal to the identity 
matrix; this reduces the -|- 3 jet production cross sec- 
tion at the LHC by less than 1%. All quarks, with the 
exception of the top quark, are considered massless. The 
top quark is considered infinitely heavy and its contribu- 
tion is neglected. The mass of the W boson is taken to be 
mw — 80.419 GeV; its couplings to fermions are obtained 
from aQED(TOz) = 1/128.802 and sin^ 6*^ = 0.230. We 
use CTEQ6L parton distribution functions for leading 
order and CTEQ6M for next-to-leading order computa- 
tions [1^, [3^. Note that we do not include the factor 



Br(H^ — > Ivi) in the results for cross-sections quoted be- 
low. 

We first discuss results for total cross sections. 
We set renormalization and factorization scales to 

= [80, 120, 160, 200, 240] GeV and calculate the cross- 
sections with the cuts defined at the beginning of this 
Section. The result of the calculation is illustrated in 
FigUI For full-color leading order cross section we find 



LO,FC 
''VK+ + >3J 



35(10) pb , 



(2) 



where the ±10 pb uncertainty from scale variation is 
shown in brackets. Calculating the same cross-section 
in the leading color approximation, we find the leading 
color adjustment parameter 

^ = C++>3,/C++>3, = 0-940(5) , (3) 

where the uncertainty indicates changes in this ra- 
tio that we observe when we change factoriza- 
tion/renormalization scales chosen in leading order com- 
putations or cuts on the final state particles. We also find 
that the TZ ratio for the W~ production is the same as 
for the . Since TZ does not depend in any significant 
way on the details of the process, applied cuts and chosen 
scales, we use the central value for TZ given in Eq.Q in 
what follows. 

At NLO we obtain the adjusted leading- color inclusive 

NLO,LC / 



NLO,aLC/. n ^ NLO,LC /• n 

cross-section, cr^^+_,:>3j(incl) = IZ ■ o•^y+_|:>3 .(mcl). 



NLO,aLC^- n 



32.4(1.5) pb. 



(4) 



This result implies (see Fig[T]) that for our choice of 
cuts and input parameters, NLO QCD corrections to 
the inclusive cross-section are very moderate for /i 
140 — 160 GeV. We also observe a remarkable reduction 
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FIG. 5: The transverse momentum distribution of the sec- 
ond hardest jet in + 3 jet production at the LHC. All 
cuts and parameters are described in the text. The leading 
color adjustment procedure is applied. The LO distribution is 
rescaled by constant factor, to ensure that the LO and NLO 
normalizations coincide. 



FIG. 6: The transverse momentum distribution of the third 
hardest jet in + 3 jet production at the LHC. All cuts 
and parameters are described in the text. The leading color 
adjustment procedure is applied. The LO distribution is 
rescaled by constant factor, to ensure that the LO and NLO 
normalizations coincide. 



in scale dependence from more than ±30% at leading 
order to only ±5% at NLO. While corrections to the ex- 
clusive cross-section are larger for similar values of /i, the 
scale independence of the exclusive NLO cross-section is 
similar to the inclusive one. In Figl^l the cross-section 
for -\- 3 jet production is shown in dependence on 
the factorization and renormalization scales. The cross 
section is smaller in this case, while the stabilization of 
scale dependence that occurs at next-to-leading order is 
very similar for W~ and production cross-sections. 

Given that NLO QCD corrections to the total cross 
sections are small, it is tempting to surmise that the 
corrections to kinematic distributions should also be in- 
significant. As we will now show, the actual situation 
is more complex. We consider kinematic distributions 
for the inclusive + 3 jet production. We choose 
to show the NLO distributions for the dynamical scale 

fj.0 = \Jp'tw~^ itt-w^ where pt,w is the transverse mo- 
mentum of the W boson as done e.g. in j33|. We note 
that for such a scale the LO cross-section is o'^y+4.>3j = 
37.6 pb and the adjusted leading color NLO cross-section 
is '^y^+'^^^j = 34.2 pb, consistent with Eq. (j4]) within 
the indicated uncertainties. The radiative corrections to 
14^ -I- 3 jet production cross-section at scale fj,o are there- 
fore small, about —10%. For the following discussion, 
scale choices in NLO computations are not very impor- 
tant since, as it turns out, shapes of NLO distributions 
are fairly insensitive to them. 

We begin by studying the transverse momentum dis- 
tribution of the leading jet. In Fig. [3] we compare NLO 
and LO predictions for scale /xq- We find that the NLO 
QCD corrections change the shape of this distribution - 
the leading order distribution underestimates the NLO 
result at small values of the transverse energy by about 
30 percent and systematically exceeds the NLO result for 



higher values of the transverse energy. A similar feature 
is observed in other distributions related to jet transverse 
momenta if the NLO result is compared to LO predictions 
with the scale /iq. 

The origin of these shape changes was recently dis- 
cussed in Ref. |38| using soft-coUinear effective theory 
and in Refs. Q i'^ connection with NLO QCD com- 
putations for -|- 3 jets production at the Tevatron and 
the LHC. Here, we recapitulate the explanation of the 
inadequacy of the scale fio given in Ref. Q. This inad- 
equacy is related to two facts: a) in the region where 
the jets have large transverse momentum, the M^-boson 
transverse momentum spectrum is softer than that of the 
jets; b) the probability of parton branching is determined 
by the relative transverse momentum of the two daugh- 
ter partons produced in that branching; such transverse 
momentum should be the appropriate scale for the strong 
coupling constant. When these two facts are combined, 
one is led to the conclusion that in the kinematic region 
where the jets have large transverse momenta, the use of 
cks(/io) in LO computations overestimates the cross sec- 
tion. At next-to- leading order, the appropriate scale for 
the strong coupling constant /i ~ ptj ^ A*o is generated 
dynamically and the cross section in that region becomes 
smaller. 

Is it possible to account for the shape modifications by 
more sophisticated LO computations? The affirmative 
answer to this question was given in Refs. [ZjiSSji where 
particular choices of scales set by e.g. the hadronic in- 
variant mass or total transverse energy in an event, were 
advocated. It should be emphasized, however, that the 
idea to employ scales of the strong coupling that are de- 
termined from local kinematics on an event-by-event ba- 
sis is not new since it is central to both parton showers 
and advanced leading order computations that employ 
matrix elements and parton shower matching (39| . 
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FIG. 7: The rapidity distribution of the leading jet for 
+3 jet inclusive production cross section at the LHC. All 
cuts and parameters are described in the text. The leading 
color adjustment procedure is applied. The LO distribution is 
rescaled by constant factor, to ensure that the LO and NLO 
normalizations coincide. 



FIG. 8: The total transverse energy distribution for W +3 jet 
inclusive production cross section at the LHC. All cuts and pa- 
rameters are described in the text. The leading color adjust- 
ment procedure is applied. The LO distribution is rescaled 
by constant factor, to ensure that the LO and NLO normal- 
izations coincide. 



Since such "local" scales capture the kinematics of 
complicated events correctly, it is conceivable that they 
produce shapes that are close to exact NLO results. We 
show the comparison of the NLO prediction with two 
leading order results in FigU) 

One LO distribution is obtained by following the MLM 
procedure whose ap phc ation to -I- 3 jet production 
is described in ReL 37]. The MLM procedure and its 
close relative the CKKW algorithm 39| are the most ad- 
vanced techniques available currently for leading order 
predictions, so it is interesting to see how it compares 
with NLO computations. We use Alpgen [i^l to gener- 
ate unweighted events that are matched to the Herwig 
[4l| parton shower. We produce hard events with up to 
five QCD partons in the final state with Alpgen, using a 
transverse momentum cut of ptj,mm = 20 (25) GeV and 
a separation parameter drj = 0.35 (0.45) To shower 
the hard events with Herwig we used i?cius = drj and 
-E't.cius = Ptj mill as matching parameters for the MLM 
prescription . We find that results are fairly indepen- 
dent of the cuts used in the generation of the hard events 
and that samples with five hard partons contribute lit- 
tle. This indicates that hard samples with yet higher 
multiplicity can be safely neglected. 

The other LO prediction shown in Fig. [4] is our im- 
plementation of the local scales in the strong coupling 
constant; it is close in spirit to the re-weighting part of 
the CKKW procedure [33]. To this end, for a given LO 
partonic event that passes jet cuts, we cluster partons 
according to the measure given by fcj_-jet algorithm^. A 



^ We note that the jet cuts can be defined with any jet algorithm; 
the fex-algorithm is only used to reconstruct the event branching 
history. 



repeated clustering gives us a "branching history" that 
can be associated with the event; at each branching the 
scale of the strong coupling constant is chosen as the rel- 
ative momentum of two daughters in the branching. We 
will refer to scales of the strong coupling constant chosen 
by this algorithm as "local" scales. Note that this pro- 
cedure is strictly a simple way to set scales of the strong 
coupling constant to reasonable values in -I- 5 parton 
leading order matrix elements. In doing so, we do not 
try to combine matrix elements of different multiplicities 
nor do we attempt to shower leading order partonic con- 
figuration. Differences between distributions produced 
with Alpgen and with the local scale procedure give an 
idea of the importance of the parton shower and Sudakov 
re- weighting. 

We point out that such modifications of leading or- 
der computations may lead to large changes in the cross- 
sections. For example, Alpgen cross-section is ^ 22 pb 
and the local scale cross-section is ~ 47 pb, to be com- 
pared with ^ 33 pb NLO cross-section. However, the 
normalization of cross-sections is a hard problem where 
next-to-leading computations or direct normalization to 
data are the only known solutions. To separate issues 
of normalization from the shape, we normalize all lead- 
ing order results in Fig. 3] to the NLO cross-section. We 
observe that both the Alpgen-fHerwig distribution and 
the local scale distribution describe the NLO result fairly 
well. Also, the proximity between the shapes of the two 
leading order results tells us that parton shower does rel- 
atively little to alter the shape of the distribution. 

We find that these observations are generic: leading 
order computations obtained with either Alpgen-|- Herwig 
or local scales are similar and they work reasonably well 
in reproducing shapes of NLO distributions. We believe 
this is important conclusion, especially in the case of 
Alpgen-|-Herwig since those programs are used by experi- 
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FIG. 9: The transverse momentum distribution of the charged 
lepton for + 3 jet inclusive production cross section at 
the LHC. All cuts and parameters are described in the text. 
The leading color adjustment procedure is applied. The LO 
distribution is rescaled by constant factor, to ensure that the 
LO and NLO normalizations coincide. 



menters as tools for understanding properties of W+ jets 
process. In order to avoid too busy plots, we choose to 
only show leading order results computed with the "lo- 
cal" scale choice for the strong coupling constant in what 
follows. We stress that for all distributions the normal- 
ization of the leading order cross-section is adjusted to 
agree with next-to- leading order result. We show the 
distribution in transverse energy of the second-hardest 
and third-hardest jet in Figs. I5I6I rapidity of the hard- 
est jet in Fig. [7| and the distribution in total transverse 
energy 7?T,tot = J2jcts j | + Pt,\ + #t in Figgl We 
also show leptonic distributions in Figs. [51 and [TUl where 
we plot the lepton transverse momentum and the miss- 
ing transverse momentum, respectively. As stated, in all 
considered cases local scales reproduce shapes of the dis- 
tributions quite well. 



III. W + 3 JET PRODUCTION AS A MODEL 
FOR BACKGROUND TO SUPERSYMMETRIC 
SEARCHES 

In this Section we investigate QCD corrections to 14^-1-3 
jet production at the LHC for a set of cuts appropriate 
in supersymmetric searches. By construction, these back- 
ground cuts seek to suppress the production of W bosons 
in association with jets as much as possible, effectively 
driving W + jet production to corners of the available 
phase-space. It is therefore unclear if QCD radiative ef- 
fects in those regions of phase-space are similar to QCD 
corrections to the production cross-sections discussed in 
the previous Section. To answer this question, we dis- 
cuss two types of cuts, very similar to those suggested 
by the ATLAS and CMS collaborations, in their planned 
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FIG. 10: The missing transverse momentum distribution for 
+ 3 jet inclusive production cross section at the LHC. All 
cuts and parameters are described in the text. The leading 
color adjustment procedure is applied. The LO distribution is 
rescaled by constant factor, to ensure that the LO and NLO 
normalizations coincide. 



searches for supersymmetry at the LHC. ^ All the input 
parameters are the same as in the previous Section ex- 
cept that we use merging parameter / = 0.7 to define 
jets using SIScone algorithm. 



A. ATLAS setup — faking jets from r decays 

We begin by considering cuts employed by the AT- 
LAS collaboration to search for SUSY with _R-parity con- 
servation. In that case, the typical signal comes from 
gluino pair-production. If each gluino decays into two 
jets and a neutralino, a SUSY signature will involve 4 
jets and missing transverse energy. A dominant back- 
ground to this process comes from Z + 4 jet produc- 
tion, with the subsequent decay of the Z-boson into 
two neutrinos. Another important background comes 
from + 3 jet production**, followed by the decays 
— > rvr — > VtI^t + hadrons, so that hadrons from 
semileptonic decay of the r lepton produce the fourth 
jet. 

One can use peculiar kinematic properties of the fourth 
jet to connect it to r decays and then reject such events 
but, because of limited efficiency in identifying r decays 
and because the cross section for -f 3 jet production is 
almost two orders of magnitude larger than the Z + 4, jet 



^ We point out that we kept cuts very similar to those used in 
the experimental studies done at 14TeV despite the fact that we 
use 10 TeV as center-of-mass energy. As a consequence cross- 
sections in this section are very small. A more realistic study 
would require adapting those cuts to the centre-of mass energy, 
but this is beyond the scope of this paper. 

^ Clearly, there is also a similar background from W~ + 3 jet pro- 
duction but we do not consider it here. 
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production cross section, it is important to consider this 
source of the background as well. 

We begin by listing a typical set of cuts that the AT- 
LAS collaboration applies to suppress the W ^ t + 3j 
background [13, [28l. l29l|. First, all jets are required to 
have transverse momenta larger than 50 GeV and the 
transverse momentum of the leading jet should exceed 
100 GeV. Second, missing energy in the event should sat- 
isfy fJr > max(100 GeV, 0.2iJT) with Ht = Y.jPt.j + 
pT- Third, no leptons with transverse momenta higher 
than 20 GeV should be present. Fourth, jets should be 
central |ryj| < 3. Finally, the event is required to be 
spherical and the cut 5't > 0.2 is applied on the trans- 
verse sphericity. We will not employ the sphericity cut in 
what follows because this observable is not coUinear safe 
at parton level. In addition, since we consider semilep- 
tonic decays of the t lepton, no high-pT lepton is present 
in our events and we do not need to employ a 20 GeV 
lepton cut. The primary observable is the distribution 
in the effective mass defined above and the range of 
a particular interest for SUSY searches, given existing 
bounds on gluino masses, is ^ 1 TeV. 

A clear exposition of the effect that the ATLAS cuts 
have on W T + 3j background at leading order was 
recently given in Ref. ^.29]. It turns out that these cuts 
primarily change the normalization of the background 
but do not significantly affect the shape of the effective 
mass distribution, especially in the region Ht ~ 1 TeV. 
We would like to understand the impact of NLO QCD 
corrections to — > r-f 3 jet on the Ht distribution. Our 
implementation of radiative corrections incorporates W 
decay to any leptonic final state but subsequent hadronic 
decays of the r-lepton are not included. Yet, as we will 
argue now, this is not necessary if all we need is an esti- 
mate of the QCD effects. 

We note that, given the above cuts and, in particular, 
the cut on the missing transverse energy, the r lepton 
produced in W decays will be highly boosted and its de- 
cay products will be very coUimated. We then completely 
neglect the angular distribution of the r decay products 
and assume a perfect collinear splitting. If, in addition, 
we neglect all the spin correlations in r decay r —^ VrQiQj , 
we conclude that the neutrino has to carry away about a 
third of the r momentum while the hadronic jet formed 
by a quark and an anti-quark from r decay has to carry 
away two-thirds of the original r momentum. We also 
expect that, since the r lepton is highly boosted, all its 
hadronic decay channels will contribute to the same jet, 
making the inclusive treatment of jet properties a rea- 
sonable approximation. 

We can implement this set up in our calculation by 
producing a W boson and letting it decay to a massless 
lepton and a massless neutrino. We then carry through 
all the steps required for the NLO QCD computation 
until the moment when the kinematics of events is ex- 
amined and weights, relevant for various histogram bins, 
are calculated. At this point, we assign one-third of the 
lepton momentum to additional missing energy carried 
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FIG. 11: Distributions in effective mass Ht for {W^ — > f)-|-3 
jet sample for ATLAS SUSY cuts described in the text. All 
cuts and parameters are described in the text. The leading 
color adjustment procedure is applied. The large difference 
between LO and NLO distributions can be absorbed by re- 
scaling the LO distribution by a constant factor. 



away by and two-thirds of the lepton momentum to 
the fourth (r) jet in the event. Note that we do not apply 
the jet algorithm to check whether or not the hadronic jet 
from T decay is sufficiently separated from the other three 
jets. ^ Since this step is not necessary for infra-red safety, 
we feel that it is entirely justified to omit it, given the 
approximate nature of our analysis. For next-to- leading 
computations, we use the leading color adjustment pro- 
cedure; we find that TZ = 0.93 is an appropriate value of 
the re-scaling parameter for ATLAS cuts. 

The results of our computation are presented in Fig. [11] 
where the LO H± distribution for our default local scale 
compared to the NLO distribution for the factorization 
and renormalization scales set to /xq. We point out that 
the shape of the leading order distribution is similar to 
that obtained with Alpgen presented in Ref. [2^, espe- 
cially at high values of the effective mass. At lower val- 
ues of the effective mass, there is a dependence on the 
modeling of r ^ hadrons transition and, given the very 
approximate nature of our procedure, it is not surprising 
that it tends to fail. It is reassuring, however, that our 
procedure seems to work quite well for high values of the 
effective mass. 

As follows from Fig. [11] the Ht mass distribution re- 
ceives large positive QCD corrections for ATLAS cuts. 
Note that distributions for local scales are not normal- 
ized to match the NLO distribution there. We studied 
the scale dependence of the leading order predictions by 
varying local scales around the central value by a factor 
of two. While we observe large ^ ±50% scale dependence 
in the LO result, the NLO QCD corrections are ~ 100% 



^ We did however impose a separation Rij = 0.5 between the r- 
lepton and the jets. 
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and are thus considerably larger than what the LO scale 
variation suggests. The scale dependence of the i/^ dis- 
tribution does decrease considerably at NLO. We find 
that NLO QCD effects provide a universal enhancement 
of Ht distribution without distorting its shape. Inter- 
estingly, the cuts on jets and missing energy presented 
at the beginning of this Section have a similar impact on 
the {W ^ t) + 3 jet background - each of the individual 
cuts reduces the magnitude of {W t) + 3 jet by a fac- 
tor between three and four, without affecting the shape 
of the Ht distribution [29|. NLO QCD effects therefore 
are comparable to the effects of the cuts and work in the 
opposite direction. 

We emphasize that, had we chosen scale fiQ also in LO 
computation, we would observe large positive NLO QCD 
effects for Hj_ distribution, in sharp contrast with large 
negative corrections for such scale choice in high-p^ ^ re- 
gions, described in the previous Section (see Fig. [3]). This 
is not surprising since, in contrast to -I- 3 jet signal 
cuts, ATLAS cuts require large amount of missing en- 
ergy, which forces W transverse momentum to be com- 
parable or larger than transverse momenta of hard jets in 
the event. Jet branching on the other hand, can occur at 
lower relative transverse momenta. Taking the relative 
transverse momentum as the correct scale for the strong 
coupling constant, it is natural that LO cross sections 
for /i — /io strongly underestimate the Ht distribution. 
This is indeed what we see when LO and NLO results 
are compared. 

We believe that this discussion shows explicitly how 
problematic extrapolation from signal to background re- 
gion can be since the NLO QCD effects for ATLAS cuts 
have no relation whatsoever to the NLO QCD effects for 
the total cross section. This mismatch happens because 
the kinematic region selected by ATLAS cuts gives negli- 
gible contribution to the total cross section. On the other 
hand, it appears that one can use low Ht < 1 TeV bins 
for ATLAS cuts to fix background normalization since 
QCD effects seem to be //T-independent and SUSY con- 
tamination in \ow-Ht bins is small. 



B. CMS indirect lepton veto cut 

How robust is the situation discussed in connection 
with ATLAS cuts? To answer this question, we study 
another example of background cuts. Those cuts are 
adopted by the CMS collaboration for SUSY searches 
at the LHC [13, [3l|. The target signal is gluino pair 
production and the final state involves jets and missing 
transverse energy. 

The CMS collaboration does not veto leptons directly. 
Rather, cuts are designed in such a way that the contri- 
bution of VF-l-jets becomes naturally small. Such cuts are 
usually referred to as indirect lepton veto cuts. We ap- 
proximate the CMS indirect lepton veto cut by requiring 
that there are three or more jets in the event. The miss- 
ing energy in the event should be large, i^miss > 200 GeV. 




600 800 1000 1200 1400 1600 
Ht,24 [GeV] 



FIG. 12: Distributions in reduced transverse mass Ht,24, for 
-I- 3 jet events with CMS SUSY cuts that define indi- 
rect lepton veto procedure as described in the text. All cuts 
and parameters are described in the text. The leading color 
adjustment procedure is applied. 



The leading jet in the event should be very central 
|??ioad jetl < 1-7 while all other jets should be in the cen- 
tral region |77othcr jctsi < 3. Jets are defined with the 
transverse momentum cut of pT,j > 30 GeV but the 
transverse momentum of the leading and sub-leading jets 
should be larger that 180 and 110 GeV, respectively. Lep- 
tons from W decays should satisfy the same cuts as jets 
but lepton transverse momentum can not be the largest 
or next- to-largest in a particular event; experimentally, 
this requirement is implemented by cutting on the frac- 
tion of electromagnetic energy carried by a "jet" . Fi- 
nally, a particular effective mass is required to be large 

4 

Htm = J2 Pt,j + -Eniiss > 500 GeV. To calculate the 

sum in this formula one orders leptons and jets accord- 
ing to their hardness, disregards the leading jet and sums 
over transverse momenta of second-to-leading, third-to- 
leading and fourth-to-leading particles/jets. 

We show the result of our computation of the NLO 
QCD corrections to the W -I- 3 jet cross section in case 
of CMS-style cuts in Figs. 1121131 where distributions in 
iJT,24 and missing energy are plotted. We again use /io, 
as the factorization and renormalization scales and vary 
it by a factor two up and down to estimate scale uncer- 
tainties. The NLO corrections for these cuts change the 
LO result by —40% to —10% depending on the scale cho- 
sen in LO computations. For fi — ^q, the corrections are 
about —10% and no significant changes of shape are ob- 
served. In this case, the scale variation at leading order 
gives a good indication of the size of NLO QCD correc- 
tions. 

It is striking that the magnitude of NLO QCD correc- 
tions for CMS cuts is in strong contrast with the magni- 
tude of NLO QCD effects for ATLAS cuts, discussed in 
the previous Section. This emphasizes the dependence of 
NLO QCD corrections on exact implementation of kine- 
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Pt, miss [GeV] 

FIG. 13: Distributions in the missing transverse energy for 
W+ + 3 jet events with CMS SUSY cuts that define indi- 
rect lepton veto procedure as described in the text. All cuts 
and parameters are described in the text. The leading color 
adjustment procedure is applied. 



matic cuts even if such cuts are designed to target very 
similar physics beyond the Standard Model. On the other 
hand, we find that shapes of basic distributions employed 
in supersymmetric searches are described fairly well by 
leading order computations, for both ATLAS and CMS 
cuts. If one can verify that, say, low-_ffj^ bins are not 
contaminated by New Physics, those bins can be used to 
determine the normalization of the background. 

IV. CONCLUSIONS 

We have discussed the NLO QCD corrections to M^ + S 
jet production at the LHC. We found that the inclusion 
of NLO QCD corrections leads to a significant reduction 
in dependence of LO results on the renormalization and 
factorization scales; the residual uncertainty associated 
with the total cross section is ±5%. We showed that 
small corrections to total cross sections do not necessar- 
ily imply that corrections to differential distributions are 
small and there is a high degree of non-uniformity in 
these corrections across the available phase-space. 

It should be stressed that the last statement depends 
upon renormalization and factorization scales chosen in 
leading order computations. In particular if leading order 

calculations are done with the scale /i = \ p^^^ + ™vi/ 



we find a large difference in shapes between LO and NLO 
distributions. On the other hand, it is clear a'priori that 
better results are achievable if scales are chosen based on 
local probabilities for jet branching. Here we have shown 
explicitly that when a local scale choice for the strong 
coupling constant is employed in leading order compu- 
tations, such computations reproduce shapes of various 
NLO distributions quite well. Note that any leading or- 
der computation matched to parton shower in the spirit 
of CKKW procedure [3§| does employ such local scales 
and our NLO analysis therefore confirms that, as far as 
shapes of various kinematic distributions are concerned, 
this is a very reasonable procedure. 

The production of M^-bosons in association with three 
jets is an important background for SUSY searches in 
jets -I- missing energy channels. We studied NLO QCD 
corrections to cuts employed by ATLAS and CMS col- 
laborations for SUSY searches and found that such cor- 
rections are not at all correlated with corrections to the 
total cross sections. It is peculiar that the magnitude of 
NLO QCD corrections to, say, effective transverse mass 
distributions, is very different for ATLAS and CMS cuts 
in spite of the fact that these cuts are designed to serve 
the same purpose. We find large (~ 100%) corrections 
for ATLAS and smaU (- 10%) QCD corrections for CMS 
cuts. We believe that this non-uniformity of corrections 
and their apparent strong dependence of the experimen- 
tal set-up emphasizes the need for extending NLO QCD 
studies to other relevant backgrounds such as -I- 4 jets 
and Z + 3,A jets. We hope that techniques for NLO 
QCD computations developed in recent years make such 
computations possible. 
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